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ABSTRACT 

The ellipticity function of cosmic voids exhibits strong dependence on the amplitude of the linear 
matter power spectrum. Analyzing the most recent void catalogs constructed by Foster and Nelson 
from the fifth data release of the Sloan Digital Sky Survey, we measure observationally the ellipticity 
function of giant galaxy voids. Then, we incorporate the redshift distortion and galaxy bias effect into 
the analytic model of the void ellipticity function and fit it to the observational result by adjusting the 
value of the power-spectrum normalization with the help of the generalized ^-minimization method. 
The best-fit normalization of the linear power spectrum is found to be cr 8 = 0.90 ± 0.04. Our result 
is higher than the WMAP c 8 -value but consistent with that from the recent work of Liu and Li who 
have constructed a new improved CMB map independently. 
Subject headings: cosmology: theory — large-scale structure of universe 



1. INTRODUCTION 

The normalization amplitude of the linear matter 
power spectrum is one of the key cosmological pa- 
rameters that are required to complete the theoreti- 
cal description of the initial conditions of the universe 
(jTegmark et al.l [20061 ). It is often quantified in terms 
of erg, the rms fluctuations of the linear density field 
within a top-hat spherical radius 8/i _1 Mpc. Various ob- 
servables have so far been used to constrain the valu e 
of erg: the cluster abundance (e.g. . iHenrv et afll2009f). 
the w eak lensing cosmic shear (e.g.. lVanW acrbck e et all 
2000), strong lensing arc statistics (e.g., IHuterer fc Mai 
2004), the cluster shapes (|Le3 |2006), and the cosmic mi- 



crow ave background radiation (CMB) te mperature map 
(e.g-. lDunklev et al.ll200l iLiu fc Lill2009l ). Yet, these ob- 
servables depend not solely on as but concurrently on the 
other key parameters such as the matter density parame- 
ter fi m , primordial non-Gaussianity parameter /nl, and 
dark energy equation of state w. Furthermore, it has 
been realized that complicated systematics involved in 
the measurement of these observables could bias strongly 
the estimates of erg. Hence, to break the parameter de- 
generacy and to diminish any systematic bias, it is very 
important to consider as many alternative probes as pos- 
sible. 

Recently, ()Park fc Led 120071 hereafter, PL07) have 
proposed the void ellipticity function as another sensi- 
tive probe of as- Noting that the shapes of voids are 
modulated by the competition between tidal distortion 
and cosmic expansion, they have analytically derived the 
void ellipticity function under the assumption that the 
dynamics of void galaxies can be well described by the 
Zel'dovich approximation just as that of the dark matter 
particles are in the linear regime. They have tested their 
model ag ainst the results from the Millennium Run sim- 
ulations (|Springel et al.ll2~005f ) , which proved the validity 
of the PL07 model. For the comparison with observa- 
tional data, however, the PL07 model has to be extended 
to incorporate the redshift distortion effect since in prat- 
ice the void ellipticities can be measured only in redshift 
space. 



Moreover, there is one condition that the success of 
the PL07 analytic model is contingent upon. Its va- 
lidity has been tested only for the case that the voids 
are found through the specific void-finding algorithm of 
iHovle fc Vogelevf (|2002l hereafter HV02). Since there is 
no unique definition of voids, the ellipticity distribution 
may well depend on the way in which voids are identi- 
fied (jColberg et all 2008). For the fair comparison with 
the PL07 model, the HV02 algorithm should be consis- 
tently used for the id entification of voi d s from observa- 
tions. Very recently, iFoster fc Nelson! ([2009, hereafter 
FN09) have constructed a catalog of 232 voids from the 
Sloan Digital Sky Survey Data Release 5 (SDSS DR5). 
Now that the voids of the FN09 catalog are identified 
using the HV02 algorithm, it must provide the most op- 
timal dataset against which the PL07 analytic model can 
be compared. 

In this Letter our goal is to constrain the value of as 
by comparing the extended PL07 analytic model of the 
redshifted void ellipticity function with the observational 
result from the FN07 catalog of SDSS voids. 

2. REDSHIFTED VOID ELLIPTICITY FUNCTION 

Let us first give a brief overview on the PL07 theoret- 
ical model. An analytic expression for the probability 
density distribution of the minor-to-major axial ratio, u, 
of a void at redshift z on the Lagrangian scale Rl was 
found by PL07 as 

dfip[fi. v\8 = S v ;a(z, R L )} 
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where a(z, Rl) represents the rms fluctuations of the lin- 
ear density field smoothed on scale Rl at redshift z, and 
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{?/, /i} (with v < fi) represent the two axial ratios of cos- 
mic voids that can be obtained from the inertia momen- 
tum tensors of the anisotropic spatial positions of void 
galaxies. The key concept of this analytic expression is 
that the two axial ratios, v and [i, are related to the 
largest and second to the largest eigenvalues, Ai and A2, 
of the tidal field smoothed on the scale Rl as 
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where S v denotes the critical density contrast of a void 
linearly extrapolated to z = 0. 
PL07 calculated the value of 6 V as the galaxy number 



density contrast as 5 V = (n vg 



n g )/n g where n vg and 



represent the number density of void galaxies and the 
mean number density of all galaxies in a given sample. 
PL09 found that S v « 0.9 on average but also noted a 
tendency that S v decreases gradually with the sizes of 
voids. The Lagrangian scale radius Rl was calculated 
as Rl = (1 + (5„) 1 / 3 ii£;/(l + z). Here Re represents the 
effective (comoving) spherical radius of a void defined 
as 4ttR e /3 — V with the void volume V. The values 
of S v and Re have to be determined from the observed 
voids that are to be used for comparison. It is worth 
mentioning here that this relation between Rl and Re 
holds good also in redshift space. 

Defining the ellipticity of a void as e = 1 — v, the 
probability density distribution of the void ellipticities 
on scale Rl at redshift z is calculated as p(e; z, Rl) = 
p(l — v\ z, Rl)- PL07 originally deri ved Equation s (Q][3]) 
for the present epoch z = 0. It was iLee fc Park! (2009) 
who extended the analytic model to higher redshifts , 
according to which the dependence of p(e\ z, Rl) on z 
and Rl comes from the dependence of a 2 on z and Rl- 

/oo 
A 2 (k)W 2 {kR L )d\nk, (4) 
-00 

where W(kR) is the top-hat filter of scale radius Rl, 
D(z) is the linear growth factor normalized as D(0) = 1, 
A 2 (fc) is the dimensionless linear matter power spec- 
trum. The functional shapes of D(z) and A 2 (k) vary 
with the background cosmology. Assuming a ACDM cos- 
mology, we use the following formula for D{z) and A 2 (k) 
(jLahav et al.lll99lL iBardeen et al.lll986h . 
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where q = /c/[fi m /i 2 Mpc _1 ] (jPeacock fc Doddsl [l99l 
and n s is a spectral index of the primordial power spec- 
trum. The void ellipticity function, /(e;z,i?i), is now 
defined as the differential number density of voids as a 



1 There are typos in Equations (2) in Lee & Park (2009), which 
are corrected here 
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Fig. 1 . — Analytic model of the rcdshifted void ellipticity function 
for the three different cases of the linear power-spectrum normal- 
ization (as = 0.7,0.8 and 0.9 as dotted, dashed and solid lines, 
respectively). 

function of e, N tv p(e; z, Rl) where N tv is the total num- 
ber of voids. To extend the above analytic expression 
for the redshifted void ellipticity distribution, we con- 
sider the simplest case where the redshift distortion ef- 
fect can be described by a single linear distortion param- 
eter, 0, which is related to the background cosmology as 
(Hamilton 1991, and references therein): 
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where b g is a linar galaxy bias factor measured in real 
space. Since we are interested in the redshift distortion 
effect on voids, we should use a b ias fa ctor of the void 
galaxies, say b vg . iBasilakos et al.l (|2007l ) have measured 
the real-space clustering of the HI galaxies and found 
that the galaxies in low-density regions are anti-biased, 
exhibiting b vg w 0.68. Now, the redshifted power spec- 
trum can be ap proximated at first order as (|Kaiserlll987l ; 
iHamiltorillQQSh 



(7) 



Ai(k)=[l + -f3 + -f3 2 A 2 (k), 



where /3 is calculated using the void galaxy bias factor 
b vg = 0.68. Replacing the real-space power spectrum 
A 2 (k) by the redshifted power spectrum A 2 (k), we fi- 
nally obtain an analytic expression for the redshifted void 
ellipticity function, f s (e; z, Rl)- Figure Q] plots the ana- 
lytic predictions of the redshifted void ellipticity function 
for the three different cases of the linear power-spectrum 
normalization (erg = 0.7,0.8 and 0.9 as dotted, dashed 
and solid, respectively). The Lagrangian void scale, red- 
shift, and total number of voids are set at the values con- 
sistent with the ones used for the void catalog (see §3). 
As it can be seen, the void ellipticity function depends 
very sensitively on the as value. 

3. COMPARISON WITH OBSERVATIONAL RESULTS 

FN09 extended the HV02 void-finding algorithm to 
improve its statistical robustness and applied it to the 
volume-limited spectroscopic sample of the galaxies from 
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Fig. 2. — Comparison of the observational results (solid dots) 
with the analytic model (solid line) with the best-fit value erg £S 
0.9. The errors include both the cosmic variance and the Poisson 
noise. The downward arrows represent f s + 2a for an ellipticity bin 
where f s < a. 

SDSS DR5 to construct a catalog of 232 voids (available 
at http://physics.ubishops.ca/sdssvoids). The volume 
limited sample has a total of 52281 galaxies in a volume of 
21310400 /i _3 Mpc 3 (C. Foster in private communication) 
. The mean galaxy number density of this sample is thus 
n g ps 2.45 x 10~ 3 in unit of /i 3 Mpc -3 . In the catalog is 
listed the redshift z, effective spherical radius Re, num- 
ber of the void galaxies N vg , three axis-lengths (a, b, c 
with a < b < c) of the best-fit ellipsoids of each SDSS 
void. The effective spherical radius Re of each void is 
related to its three axis-lengths as R E = abc. FN09 de- 
termined the best-fit ellipsoid of each void with the hel p 
of the prescription of Uang-Condell &: Hernquistl (|2001h . 
assuming a flat ACDM cosmology with f2 TO = 0.28, 
Six = 0.72 and Ho — 100ft.. Throughout this paper, we 
also adopt the same cosmology, setting the other key pa- 
ram eters at n s = 0.96 an d h = 0.71 from the WMAP pri- 
ors (jDunklev et al.ll2009f) . The mean values of z, Re, and 
N vg averaged over all voids are found to be z = 0.114, 
R E = 24.645 /i _1 Mpc, and N vg = 7, respectively. 

Using information on a, b, c, we measure the ellipticity 
of each void as e = 1 — c/a. Binning the ellipticities as 
£j, counting the number of SDSS voids in each ellipticity 
bin, and dividing the void number counts by efej, we de- 
termine the observational void ellipticity function, f s (ei). 
To estimate the statistical errors in the measurement of 
f s (si), we separate the voids into six subsamples each 
of which has approximately the same number of voids, 
and then calculate the void ellipticity function separately 
from each subsample. The errors at each ellipticity bin 
£j is now calculated as the standard deviation of f(£i) 
between the six subsamples, which include both the cos- 
mic variance and the Poisson noise. Figure [2] plots the 
resulting void ellipticity function as dots. 

The density contrast of each void is computed as 5 V = 
(n vg — n g )/h g where n vg = N vg /V with V = (Awabc/S). 
And the mean density contrast averaged over all voids 
is determined to be 8 V ~ —0.97. Using the values of z, 
Re and S v , the Lagrangian void scale is obtained to be 
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Fig. 3. — Contours of the reduced x 2 (Xr = b 2, 3 as thick solid, 
thin solid and dotted lines, respectively) in the f! m — <rg plane. 

R L ~ 8.534/iT 1 Mpc. We fit the analytic model of the 
rcdshifted void ellipticity function to the observational 
results by adjusting the value of a% with the help of the 
generalized \ 2 minimization method where \ 2 lS given 
as 

X 2 = [JO {£i) ^ f T [£i . as)]C ^ [f O {£j) _ f T {£i . a8)l (8) 

where ff(£i) and /J(£i; erg) denote the observed and the 
theoretical ellipticity function at the i-th ellipticity bin, 
respectively. The covariance matrix (Cy ) is defined as 
dj = (Af°(ei)Af° (ej)) where the ensemble average is 
taken over the six subsamples. Since the number of the 
subsamples (six) is larger than the number of the ellip- 
ticity bins (five), it is guaranteed that (C^) is invertible 
(|Hartlap et al.ll2007D . The best-fit a s that minimizes X 
is found to be erg = —0.897 ± 0.037 where the errors 

— 1 /2 

are calculated as \f2 (c? 2 x 2 /c?ct|) . Figure [2] plots the 
analytic model with best-fit erg, demonstrating that the 
analytic model agrees well with the observational results. 

Although a fixed value, Sl m = 0.28, is assumed for 
the construction of the FN09 void catalog, it is worth 
examining the parameter degeneracy between erg and fi m 
when the void ellipticity function is used as a probe of 
erg. Varying the values of erg and f2 m simultaneously, 
we repeat the whole fitting process. Figure [3] plots the 
contours of \r — 1; 2, 3 (thick solid, thin solid and 
dotted lines, respectively) in the Q, m — erg plane, where 
X 2 is the reduced x 2 - As it can be seen, the best-fit er 8 
decreases very mildly as f2 m increases. 

4. DISCUSSION AND CONCLUSION 

We have constrained the matter power spectrum nor- 
malization as er 8 = 0.90 ± 0.04 by comparing the void 
ellipticity function from SDSS DR5 with the extended 
PL07 model. It is intriguing to note that our result is 
higher than the WMAP5 value (er 8 = 0.796 ±0 .036) but 
consis tent with that from the recent work of iLiu fc Lil 
(2009) who casted a doubt on the accuracy of the cosmo- 
logical parameters estimated by the WMAP team. Using 
a new improved CMB map constructed by employing 
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an independent software scheme, iLiu &; L i ( 2009) have 
found a s = 0.921 ± 0.036. 

The advantage of using the void ellipticity function as 
a probe of the power spectrum normalization is that it 
is purely analytical, free from any nuisance parameters 
and ad-hoc assumptions. Besides, as shown recently by 
lLam et all (|2009f) . the void ellipticity function does not 
depend on the primordial non-Gaussianity parameter, 
unlike the other prominent probe of as, the cluster mass 
function. It depends most sensitively on as among the 
key cosmological parameters and thus it is in principle 
one of the most optimal probes of the power spectrum 
normalization. Yet, it is worth noting that there is one 
weak point about using the void ellipticity function as a 
cosmological probe. The number of observable voids at 
different rcdshifts is relatively small compared with that 
of observable clusters, which means that it tends to suffer 
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